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This  work  presents  a  novel  heat-integrated  fuel  cell  stack  system  with  methanol  reforming.  Its  configu¬ 
ration  is  composed  of  fuel  processing  units  (FPUs),  proton  exchange  membrane  (PEM)  fuel  cell  stack,  and 
heat  exchangers  (HEXs).  Well  mixed  methanol  and  oxygen  flows  in  contact  with  countercurrent  flowing 
water  dominates  the  production  of  hydrogen  at  the  exit  of  FPUs  and  influences  the  stack  temperature. 
The  heat  exchange  connections  can  enhance  the  utilization  of  energy  of  FPUs.  To  ensure  the  stable  steady- 
state  operation,  the  model-free  fuzzy  incremental  control  scheme  within  the  multi-loop  feedback  control 
framework  is  developed.  Finally,  the  proposed  system  integration  and  control  configuration  are  verified 
by  closed-loop  simulations. 

Crown  Copyright  ©  2009  Published  by  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  widely  recognized  as  one  of  the  most  promis¬ 
ing  clean  energy  technologies  of  the  future.  The  proton  exchange 
membrane  (PEM)  fuel  cell  is  quite  suitable  for  residential  or  auto¬ 
motive  applications  [1  ],  for  the  following  reasons,  (i)  It  can  operate 
at  a  relatively  low  temperature;  (ii)  it  has  relatively  high  power 
density;  (iii)  its  maintenance  is  simple.  However,  the  efficiency 
of  the  PEMFC  is  strongly  affected  by  the  unsteady  hydrogen  feed 
flow,  cell  temperature,  membrane  dehydration  and  fast  variations 
in  load.  Developing  a  source  of  pure  hydrogen  is  a  commercial  chal¬ 
lenge.  Liquid  hydrocarbons  such  as  methanol,  ethanol,  and  gasoline 
are  usually  treated  as  alternative  hydrogen-rich  fuel  streams  [2]. 
Notably,  the  methanol-fuelled  reformer  has  been  widely  studied 
for  fuel  cell  applications. 

Toward  the  development  of  a  fuel  reforming  system,  Lindstrom 
and  Pettersson  [3]  presented  a  compact  methanol  reformer  for  fuel 
cell  vehicles;  Lattner  and  Harold  [4]  developed  a  kinetic  reactor 
model  for  autothermal  reforming  of  methanol;  Stamps  and  Gatzke 

[5]  modeled  the  packed-bed  methanol  reformer;  Choi  and  Stenger 

[6]  proposed  an  integrated  methanol  reformer  system  to  study  the 
hydrogen  yield  and  economic  profit,  and  Wang  and  Wang  [7]  used 
a  thermodynamic  and  exergetic  analysis  of  a  PEMFC  stack  system 
with  methanol  reforming.  Recently,  Xu  et  al.  [8]  explored  the  energy 
efficiency  of  a  methane  reforming  system  with  heat  integration.  In 
fact,  the  methane-fuelled  reforming  system  has  been  inadequately 
incorporated  into  the  PEMFC  system  because  of  restrictions  on 
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operating  temperature  and  pollutant  emissions.  In  our  opinion,  an 
integrated  power  generation  system  should  depend  on  the  synergy 
between  heat  integration  and  control  design. 

With  respect  to  the  fuel  cell  control  problem,  Lauzze  and 
Chmielewski  [9]  provided  a  set  of  feedback  structures,  includ¬ 
ing  power,  temperature  and  relative  humidity  controllers;  Woo 
and  Benziger  [10]  showed  that  the  proportional-integral-derivative 
(PID)  control  could  accelerate  the  response  of  a  PEMFC  system  to 
satisfy  the  load  demand;  Methekar  et  al.  [11]  adopted  the  multi¬ 
input  multi-output  (MIMO)  PI  control  framework  to  control  of  the 
power  density  and  temperature  of  a  distributed  parameter  model 
of  fuel  cell  systems,  and  Wu  et  al.  [  1 2  ]  indicated  that  the  MIMO  con¬ 
trol  structure  can  ensure  the  highly  efficient  control  performance 
when  a  low-order  and  simplistic  model  of  the  PEMFC  systems  is 
considered.  However,  the  models-based  predictive  control  schemes 
are  adopted  to  improve  the  control  performance  by  using  the 
fuzzy  Hammerstein  model  [13]  or  neural  network  techniques 
[14].  Moreover,  fuzzy  control  schemes  [15,16]  and  fuzzy  neural 
networks  [17]  have  been  successfully  applied  to  many  fuel  cell 
systems. 

This  article  introduces  the  kinetics  and  modeling  of  fuel  pro¬ 
cessing  units  (FPUs),  including  a  methanol  reformer  (MR),  heat 
exchangers  (HEX),  a  water  gas  shift  (WGS)  reactor,  and  a  preferen¬ 
tial  oxidation  (PROX)  reactor.  In  the  proposed  system  configuration, 
the  inlet  methanol  flow  can  influence  the  amount  of  hydrogen 
produced,  and  the  countercurrent  flow  of  water  by  virtue  of  heat 
exchange  connections  can  improve  energy  utilization  and  reg¬ 
ulate  the  stack  temperature.  Furthermore,  the  multi-loop  fuzzy 
incremental  control  framework  ensures  the  satisfactory  tracking 
performance  and  reliability  of  the  proposed  heat-integrated  system 
according  to  the  closed-loop  simulation. 
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2.  PEMFC  stack  system 

A  Ballard  5  kW  PEMFC  stack  system  was  constructed  using  the 
empirical-type  model  [18,19],  which  consists  of  35  cells  in  series. 
Regarding  the  individual  PEMFC  stack  system  shown  in  Fig.  1 ,  pure 
hydrogen  is  fed  to  the  anode,  excess  hydrogen  gas  is  re-circulated, 
and  the  stack  temperature  is  changed  by  a  heat  exchange  system 
with  circulating  water. 

2.1.  Anode  and  cathode  flow  models 


Based  on  the  assumption  that  all  gases  are  ideal,  the  principles 
of  mole  conservation  can  be  used  to  model  the  reactant  flows  at  the 
anode  and  the  cathode: 


Va_  dPn2  =  .FC 
RT  d t  H2’l'n 

Vc  dPp2  _  FC 
RT  d t  °2 >in 


ka{Pu2  Pamb ) 
kc(P<D2  ~  Pamb) 


nJ 
2 F 
nt 
4 F 


(1) 


where  h{^  and  are  the  hydrogen  and  oxygen  inlet  mole 
flow-rates,  respectively.  Notably,  the  inlet  hydrogen  flow-rate  usu¬ 
ally  varies  with  the  amount  of  hydrogen  supplied,  and  the  oxygen 
flow-rate  can  be  determined  by  multiplying  the  air  inlet  flow-rate 
by  the  mole  fraction  of  oxygen  in  air  (21%). 


2.2.  Thermal  model 


The  energy  balance  of  the  stack  system  is  determined  by  the 
total  power  converted  by  the  fuel  H2  into  electricity,  the  power  con¬ 
sumed  by  the  electrical  load,  the  rate  of  heat  removal  by  the  coolant, 
Qcoob  and  the  rate  of  heat  loss  at  the  stack  surface,  Qioss.  The  dynamic 
model  is  described  by  a  couple  of  first-order  differential  equations: 


Tc  in  and  Tc,out  are  the  inlet  and  outlet  water  temperatures,  respec¬ 
tively.  mew  is  the  mass  flow- rate  of  water;  pw  is  the  density  of  water; 
Vw  is  the  volume  of  the  heat  exchange  system,  and  T  is  the  stack 
temperature. 

2.3.  Polarization  curve  model 

The  polarization  curve  is  generally  used  to  specify  the  relation 
between  the  cell  voltage  Vfc  and  current  density  I.  When  a  cell  deliv¬ 
ers  power  to  the  load,  the  load  voltage  E  is  reduced  by  the  voltage 
drop,  comprising  of  the  activation  overvoltage  Vact  and  the  ohmic 
overvoltage  Vohm.  The  output  voltage  of  a  single  fuel  cell  is  given  by 

Vfc  =  E-Vact-Vohm  (4) 

(i)  The  open-circuit  cell  potential  E  is,  as  determine  by  the  Nernst 
equation,  given  by 

RT 

E  =  1.229  -  8.5  x  1(T4(T  -  298.15)  +  —  ln[PH2(Po2)°'5]  (5) 

(ii)  With  respect  to  the  dynamics  of  the  activation  overvoltage  Vacu 
the  first-order  dynamic  accounting  for  the  effects  of  double 
layer  capacitance  charging  at  the  electrode-electrolyte  inter¬ 
faces  is  described  by 

d  Vact  _  J_ _ Eact  .  . 

dt  Cd/  RactCdi 

where  the  activation  resistance  Ract  =  Vact  IP  the  activation  drop 
Eact  is  defined  by 

Pact  =  Pi  +  &T  +  p3T\n{C02)  +  p4T\n(I),  (7) 

and  the  parametric  coefficients  /?i, . . .,  are  expressed  as 
=  -0.948 

=  0.00286 +  0.0002 In (Apc)  + 4.3  x  10-5ln(CHJ 
ft  =  7.6xl0-5  (8) 

ft  =  -1.93  x  10“4 

with 

Cb2  =  1.97  x1CT7Po2  exp 

(9) 

Ch2  =  9.174  x  1(T7Ph2  exp  (^7) 

Notably,  C02  is  the  oxygen  concentration  at  the  cath¬ 
ode/membrane  interface  and  Ch2  is  the  hydrogen  concentra¬ 
tion  at  the  anode/membrane  interface. 

(iii)  The  ohmic  overvoltage  is  given  by, 

Vohm  —  ^int  (10) 


dr 

Ct  =  Etot  ~  Pelec  ~  Qioss  ~  Qcool 
PwVwCpW — -q: —  =  rhcwCpw^Tcjn  —  Pc, out) 

+UA^T-Tcin+2Tc°Ut 

where 

Et0[  AH 
2  E 

Pelec  —  VtfaapI 

+  _  Q(T  —  Tamb) 

'Moss  — 

Qcool  —  ( hcond  +  hconvl) 


C T-Tc,in)-(T-Tc,out ) 

ln((T-TCffn)/(T-TC|oUt)) 


(2) 


(3) 


where  the  internal  resistance  Rint,  obtained  by  the  empirical 
analysis,  is  written  as, 

Riot  =  ^  (ID 

Afc 

and  the  membrane  resistivity  rM  is  given  by 

181 .6[1  +  0.03(//ftc)  +  0. 062(77303  )2(JMfc)2'5]  ,  s 

rM  =  [11 .866  -  3(//ftc)]  exp[4.18((T  -  303 )/T)]  (12) 

According  to  Eqs.  ( 1 )— (12),  a  single  fuel  cell  is  affected  by 
current  density,  cell  temperature,  hydrogen  and  the  partial 
pressure  of  oxygen.  Table  1  presents  parameter  definitions  and 
values  of  a  Ballard  5  kW  PEMFC  stack  system.  Since  all  cells  are 
in  series,  the  total  voltage  for  the  stack  is  given  by 


Vstack  —  35Vjc 


(13) 


922 


W.  Wu,  C.-C.  Pai  /  Journal  of  Power  Sources  194  (2009)  920-930 


Table  1 

Parameter  definitions  and  values  of  the  Ballard  MK5-E  system. 


Parameter 

Description 

Value 

Anode  volume  (m3) 

0.005 

ka 

Anode  flow  constant  (mol  s_1  atm-1 ) 

0.065 

Hydrogen  inlet  flow- rate  (mol  l-1 ) 

0.8 

Pamb 

Ambient  pressure  (atm) 

1 

Vc 

Cathode  volume  (m3) 

0.01 

kc 

Cathode  flow  constant  (mol  s_1  atm-1 ) 

0.065 

<2,m 

Oxygen  inlet  flow- rate  (mol  s-1 ) 

2 

F 

Faraday  constant  (C  mol-1 ) 

96,485 

AH 

Hydrogen  enthalpy  of  combustion 
(kj  mol-1  j 

285.5 

ct 

Thermal  capacitance  (kj  °C_1 ) 

17.9 

hcond 

Parameter  for  conduction  property  of 
heat  exchanger  (W°C_1 ) 

35.55 

hconv 

Parameter  for  convection  property  of 
heat  exchanger  (W°C_1  A-1 ) 

0.025 

Tamb 

Ambient  temperature  (°C) 

25 

Pc, in 

Inlet  water  temperature  (°C) 

25 

UA 

Stack  heat  transfer  coefficient  (WK_1 ) 

241 

Cpw 

Heat  capacity  of  water  (kj  kg-1  K-1 ) 

4.184 

Pw 

Water  density  (kg  m-3 ) 

1000 

Vw 

Volume  of  the  cooling  system  (m3) 

2.5  x  10-3 

T 

Time  constant 

2.06 

Afc 

Effective  cell  area  (cm2) 

232 

lm 

Membrane  thickness  (cm) 

178  x 10“4 

R 

Universal  gas  constant  (J  mol-1 1<-1 ) 

8.314 

Ca 

Double  layer  capacitance  (F) 

0.035  x  232 

X 

Membrane  resistivity  parameter 

12.5 

2.4.  Simulation  and  analysis 

The  Simulink™  was  used  to  simulate  the  PEMFC  stack  system 
whose  structure  is  given  in  Appendix  A.l.  Fig.  2  shows  that,  when 
the  current  is  step  changed  from  0.5  to  45  A  at  the  time  20  min, 
the  stack  temperature  obviously  arises  and  may  cause  the  mem¬ 
brane  dehydration  when  the  temperature  exceeds  170  °C.  If  ±5% 
step  changes  of  mew  are  considered  in  Eq.  (2),  then  the  variation 
of  the  mass  flow-rate  of  water  induces  the  asymmetric  responses 
of  the  stack  temperature,  as  shown  in  Fig.  3.  In  both  figures,  the 
circulating  water  flow  clearly  affects  the  temperature  of  the  stack 
system. 


Fig.  2.  Open-loop  responses  of  voltage  and  temperature  of  the  PEMFC  stack  system 
when  the  step  change  of  the  current  appears. 
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Fig.  3.  Open-loop  responses  of  the  PEMFC  stack  system:  (A)  the  step  change  of  water 
flow- rate;  (B)  corresponding  stack  temperature. 


3.  Fuel  processing  units 

3.1.  Unit  connections 


FPUs  mainly  consist  of  HEXs,  an  MR  system,  a  WGS  reactor,  and  a 
PROX  reactor.  Their  dynamics  are  simplified  using  ideal-gas  models. 
In  an  MR  system,  water  and  methanol  flows  are  well  mixed  to  pro¬ 
ceed  the  steam  reforming,  and  an  excess  amount  O2  (air)  is  added 
to  satisfy  the  demand  for  partial  oxidation  of  methanol  (POM).  The 
equation  of  continuity  for  species  in  a  multicomponent  reacting 
mixture  can  be  used  to  determine  the  exit  mole  flow-rates  of  com¬ 
ponents  nj^,  and  nj^0  according  to  specified  reaction  rates 
in  Eqs.  (B1)-(B3)  under  prescribed  inlet  conditions  and  catalyst 
weights  in  Table  2.  Since  CO  is  produced  by  direct  decomposition 
of  methanol,  a  mixture  of  CO  and  H20  takes  part  in  the  reversible 
exothermic  reaction,  known  as  a  water-gas  shift  reaction,  to  pro¬ 
duce  hydrogen  and  reduce  the  concentration  of  CO.  The  reaction  of 
the  WGS  system  is  described  in  Appendix  B  and  the  exit  mole  flow- 
rates  of  components  hj^GS,  h^GS,  h^GS  and  nj^0  are  determined  by 
the  following  material  balances: 


hWGS 

H2 

^ WGS 
C02 


CO 

„WGS 

H20 


=  "h*  +  rwes 

=  rWGS 
=  «co  -  rWCS 

=  <0  —  rwcs 


(14) 


where  the  rate  of  the  WGS  reaction  rWGS  is  determined  by  solving  Eq. 
(B4)  at  the  prescribed  operating  temperature.  The  outlet  flows  of  the 
WGS  reactor  include  a  small  amount  of  CO,  which  is  typically  less 
than  1%.  The  selective  oxidation  reaction  of  CO  in  a  PROX  reactor  can 


Table  2 

Parameter  values  for  fuel  processing  units. 


FPUs 

Catalyst  (g  cat) 

Volume  (1) 

MR 

Cu/Zn0/Al203:  100 g  (Steam  reforming) 

Cu/Zn0/Al203:  10 g  (Partial  oxidation  of  methanol) 
Cu/Zn0/Al203i  140 g  (Methanol  decomposition) 

10.2 

PROX 

Pt-Fe/7-alumina:  40  g 

10 

WGS 

Cu/Zn0/Al203:  100  g 

0.5 
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FPUs 


Fig.  4.  PEMFC  stack  system  with  FPUs  and  circulating  water  system. 


prevent  poisoning  of  the  fuel  cell  electrodes.  Appendix  B  presents 
the  kinetics  of  the  PROX  reactor. 


3.2.  System  integration 

In  the  proposed  configuration  in  Fig.  4,  the  excess  hydrogen  gas 
is  circularly  fed  into  the  PEMFC  stack;  methanol  is  usually  selected 
as  an  alternative  hydrogen-rich  fuel  stream,  and  the  countercurrent 
water  flow  design  is  employed.  In  the  approach,  the  original  circu¬ 
lating  device  for  cooling  the  stack  system  is  replaced  with  external 
circulating  water  flow.  To  ensure  the  heat-integrated  system  is  com¬ 
pact,  and  to  reduce  its  heat  duty,  the  FPUs  are  characterized  by  the 
following. 


•  If  the  redundant  ***perheater  or  post-combustion  unit  is  omitted, 
then  the  heat  exchange  connections  for  HEXs  aim  to  improve  the 
energy  utilization  against  the  heat  duty. 

•  The  countercurrent  water  flow  from  the  stack  to  the  FPUs  plays 
the  role  of  reactant  and  heat  exchanges. 

•  A  well  mixed  lowing  CH3OH  and  02  flows  comes  into  contact 
with  the  countercurrent  water  flow,  dominating  the  heat  duty 
and  production  of  hydrogen  by  each  reactor. 

•  At  the  exit  of  the  PROX  reactor,  the  amount  of  hydrogen  produced 
as  well  as  the  CO  concentration  can  satisfy  the  specification  of  the 
inlet  flow  of  the  PEMFC  system. 

The  calculation  of  the  inlet  and  outlet  temperatures  of  each  reac¬ 
tor  depends  on  the  energy  balance.  Referring  the  unit  connections 
in  Fig.  5,  the  energy  balance  of  the  MR  system  is  expressed  as 

Qf(Tf)  +  Q'"(T«)  -  -  Q20Ut(T20Ut)  +  EMR  =  0  (15) 

where  notation  T  stands  for  the  red  flow  and  notation  ‘2’  repre¬ 
sents  the  blue  flow.  EMr  represents  the  rate  of  generation  of  heat  by 
chemical  reactions  of  the  MR  system.  Assume  that  all  reacting  units 
are  denoted  as  adiabatic  steady-state  reactor  systems.  In  the  PROX 
system,  the  flowing  water  is  heated  from  298  K  at  the  inlet  to  an 
outlet  temperature  over  373  K,  and  the  latent  energy  for  gasifying 
the  water  is  given  by 

Qh"o  =  "h20X[cp2°(,)(373  -  298)  +  44010  +  C^°(s)(T°ut  _  373)1 

(16) 

Appendix  A  presents  a  Simulink™  model  of  the  heat-integrated 
FPUs.  Fig.  6  shows  that  the  inlet  and  outlet  temperatures  and  the 
amount  of  hydrogen  produced  by  each  reactor  can  be  individually 
determined  using  kinetics  and  by  modeling  of  every  unit,  as  shown 


CHpiU- - >j^  hex!  Mixer 

- ^PPi - HEX4  I 


-rf]  rCT 

j_  T  T2  _ j  HEX3 


<5  H2° 


Fig.  5.  CH3OH  and  02  flows  come  into  contact  with  the  countercurrent  water  flow  in  FPUs. 
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Fig.  6.  (A)  Specific  amount  of  hydrogen  produced  at  each  reactor  of  FPUs;  (B)  corresponding  temperature  at  each  reactor  of  FPUs. 
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Table  3 

Parameter  values  for  heat  exchanger. 


in  Appendix  B.  In  Fig.  6(A),  1  mole  CH3OH  with  1.2  mole  02  (red 
line)  is  fed  into  the  mixer  point  and  2  mole  H20  as  the  coolant  and 
reactant  (blue  line)  is  simultaneously  fed  into  the  PROX.  The  pro¬ 
duced  hydrogen  (over  2.5  mole)  exits  the  PROX  reactor  as  the  inlet 
fuel  of  the  PEMFC  stack  system.  In  Fig.  6(B),  the  inlet  and  outlet  tem¬ 
peratures  of  FPUs  are  evaluated  with  reference  to  specifications  of 
heat  exchangers  and  the  thermodynamic  properties  of  the  reactants 
listed  in  Tables  2-4.  Notably,  Eqs.  (B7)-(B9)  describe  the  dynamics 
of  FIEX  for  multicomponent  inlet  and  outlet  flows.  However,  the 
outlet  temperature  of  HEX4  exceeds  80  °C  which  is  inadequate  for 
the  feed  flow  of  the  PEMFC  stack  system.  Therefore,  an  extra  coolant 
device  is  used  to  cool  rapidly  the  hydrogen  flow. 

4.  Control  implementation 

Although  flowing  methanol,  oxygen  and  water  are  reactants  of 
MR,  WGS  and  PROX  reactors,  Fig.  7(A)  shows  that  step  changes  of 
the  inlet  water  flow-rate  obviously  change  the  stack  temperature 
and  Fig.  7(B)  indicates  that  step  changes  in  the  inlet  methanol  flow- 
rate  strongly  change  the  hydrogen  pressure  at  the  anode.  Notably, 
both  flows  can  be  treated  as  manipulated  inputs. 

4.1.  Fuzzy  incremental  control  design 

The  model  of  the  proposed  heat-integrated  power  generation 
system  is  a  simplistic  combination,  because  the  dynamic  behav¬ 
ior  of  FPUs  is  ignored  and  the  PEMFC  is  an  empirical-type  model.  In 
fact,  the  modules  in  Appendix  A  are  incomplete  because  Simulink™ 
blocks  do  not  consider  the  effects  of  the  condition  of  the  catalyst, 
mass  and  heat  transfers  at  the  electrodes,  the  operational  limits 
of  the  power  plant,  and  balance  of  plant  (BOP)  on  the  electro¬ 
chemical  efficiency  [22]  are  not  involved  in.  Additionally,  problems 
of  unmodeled  dynamics  and  model  errors  are  inevitable.  In  our 
approach,  a  model-free  fuzzy  incremental  control  is  employed 
because  of  the  PID-like  configuration  and  the  sensible  tuning  gains. 
Fig.  8(A)  depicts  the  fuzzy  incremental  control  blocks,  in  which  the 
PD  controller  output  un  is  expressed  as 

ui=KpEi  +  KdCEi  (17) 

where  I<p  and  I<d  are  tuning  parameters,  and  error  Fz  and  change  in 
error  CFZ  represent  the  inputs  of  fuzzy  logic  controller.  Using  the 
backward  difference,  the  change  in  error  is  defined  by 

CEj  =  Ei  ~ T'"1  (18) 

Is 

where  Ts  is  the  sampling  period.  The  fuzzy  logic  controller  consists 
of  the  fuzzy  rule  (Table  5)  and  the  A -type  membership  func- 


Table  4 

Heat  capacity  of  components  of  fuel  processing  units. 


Component 

Heat  capacity  (J  kg-1  K-1 ) 

CH3OH 

61.43 

02 

29.4 

H20(d 

75.3 

H20(g) 

33.6 

CO 

29 

h2 

28.8 

co2 

43.81 

Fig.  7.  Open-loop  responses  of  the  heat-integrated  power  generation  system  with 
respect  to  step  changes  of  CH3OH  and  H20  flows:  (A)  stack  temperature;  (B)  hydro¬ 
gen  pressure  at  the  anode. 

tion/A.  The  control  structure  that  is  shown  in  Fig.  8(B)  has  been 
proven  to  generate  good  results  for  control  applications  as  well 
as  to  be  easily  implemented  in  hardware.  When  two  inputs  and 
one  output  membership  were  scaled  in  the  [-1+1]  interval,  the 
input-output  mappings  is  specified  by  a  fuzzy  inference  mecha¬ 
nism  with  respect  to  3  x  3,  5  x  5  and  7x7  linguistic  variables  as 
displayed  in  Fig.  9(A-C),  respectively.  The  corresponding  control 
surfaces  are  more  or  less  bumpy.  In  Table  5,  the  fuzzy  set  of  linguis¬ 
tic  variables  is  defined  as  follows;  NB  (Negative  Big),  NM  (Negative 
Medium),  NS  (Negative  Small),  ZE  (Zero  Error),  PS  (Positive  Small), 
PM  (Positive  Medium)  and  PB  (Positive  Big).  Table  5  also  indicates 
the  IF-THEN  relationship  between  the  membership  functions  with 
respect  to  the  linguistic  variables  for  both  inputs  Et  and  CFZ.  More- 


Table  5 

Symmetric  arrangements  of  fuzzy  rule  base  with  respect  to  3x3 
(blue),  5x5  (red  and  blue)  and  7x7  linguistic  variables. 
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Fig.  8.  (A)  Simulink™  blocks  for  the  fuzzy  incremental  controller;  (B)  fuzzy  logic  control  structure. 


over,  the  output  uf  of  the  fuzzy  logic  controller  is  written  as 

V{  =fA(KpEi,I<dCEi )  (19) 

and  the  fuzzy  incremental  controller  output  is  expressed  as 

n 

U„  =  Kw  *  YJa(Kp£„  I<dCEi)*  Ts  (20) 

2=1 


where  I<w  is  the  proportional  gain.  Based  on  the  developments  of 
PID-fuzzy  control  techniques  [23-25],  the  above  structure  is  similar 
to  a  fuzzified  proportional-integral-derivative  (PID)  controller. 

4.2.  Control  effectiveness 

To  evaluate  the  effectiveness  of  PID-fuzzy  control,  a  single-input 
single-output  (SISO)  control  framework  based  on  the  manipulation 


Fig.  9.  Plots  for  the  input-output  mappings  with  respect  to  fuzzy  sets:  (A)  3  x  3;  (B)  5  x  5;  (C)  7x7. 
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Fig.  10.  SISO  fuzzy  control  architecture  by  tuning  parameter  I<w:  (A)  stack  temper¬ 
ature;  (B)  hydrogen  pressure  at  the  anode;  (C)  corresponding  manipulation  of  the 
water  flow-rate. 


(C)  time(min) 


Fig.  11.  SISO  fuzzy  control  architecture  by  alternative  fuzzy  rule  bases:  (A)  stack 
temperature;  (B)  corresponding  hydrogen  pressure  at  the  anode;  (C)  corresponding 
manipulation  of  the  water  flow-rate. 


c 

(D 

CD 

O 


(B) 


Fig.  12.  SISO  fuzzy  control  architecture  by  tuning  parameter  I<w  for  the  perturbation 
of  methanol  flow  at  the  inlet:  (A)  stack  temperature;  (B)  hydrogen  pressure  at  the 
anode;  (C)  corresponding  manipulation  of  the  water  flow-rate. 


Fig.  13.  SISO  fuzzy  control  architecture  by  alternative  fuzzy  rule  bases  for  the  pertur¬ 
bation  of  methanol  flow  at  the  inlet:  (A)  stack  temperature;  (B)  hydrogen  pressure 
at  the  anode;  (C)  corresponding  manipulation  of  the  water  flow-rate. 
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of  the  single  inlet  water  flow-rate  is  considered.  For  fixed  parame¬ 
ters  ( Kp ,  I<d)  =  (0.1, 30),  Fig.  10(A)  shows  that  when  the  load  demand 
is  initially  changed  from  40  to  60  A,  the  stack  temperature  can  be 
asymptotically  regulated  at  75  °C  by  tuning  I<w,  but  the  hydrogen 
pressure,  plotted  in  Fig.  10(B),  clearly  decrease.  Similarly,  Fig.  11(A) 
shows  that  the  stack  temperature  can  be  asymptotically  regulated 
at  75  °C  using  alternative  fuzzy  rule  bases.  Notably,  the  fuzzy  sets 
5x5  and  7x7  ensure  superior  output  regulation,  but  the  high-gain 
(rapid)  control  responses,  shown  in  Fig.  11(C),  appear.  Fig.  12(A) 
shows  that  when  the  step  change  of  inlet  flow-rate  of  methanol, 
from  1  to  3mols-1,  occurs,  the  stack  temperature  can  be  main¬ 
tained  at  the  desired  level  by  tuning  I<w,  but  the  hydrogen  pressure, 
shown  in  Fig.  12(B),  clearly  increases.  Moreover,  Fig.  13  shows  that 
the  output  regulation  of  the  stack  temperature  is  greatly  improved 
by  the  use  of  alternative  fuzzy  rule  bases.  According  to  the  afore¬ 
mentioned  closed-loop  simulations,  the  stack  temperature  can  be 
regulated  by  manipulating  the  water  flow-rate,  but  the  hydrogen 
pressure  at  the  anode  is  very  sensitive  to  load  changes.  Notably,  the 
countercurrent  water  flow  fails  to  dominate  the  rate  of  hydrogen 
production.  The  SISO  fuzzy  control  scheme  is  characterized  by  the 
following. 

•  The  alternative  fuzzy  rule  bases  can  improve  the  output  regu¬ 
lation  performance,  but  they  clearly  increase  the  computational 
burden. 


•  A  single  controller  parameter  by  tuning  I<w  is  reliable. 

•  The  SISO  control  configuration  fails  to  solve  the  problem  since  the 
rate  of  hydrogen  production  clearly  declines. 

To  overcome  the  drawback  of  the  SISO  control  framework,  two 
fuzzy  incremental  controllers  are  implemented  to  regulate  the  stack 
temperature  and  the  hydrogen  pressure  of  the  anode  by  manipu¬ 
lating  the  inlet  flow-rates  of  both  water  and  methanol.  Based  on 
the  assumption  that  the  control-loop  interaction  can  be  decoupled 
by  the  selection  of  input-output  pairings,  a  multi-loop  SISO  control 
configuration  is  considered.  According  to  the  fuzzy  control  struc¬ 
ture  in  Eq.  (20),  the  controller  output  for  inlet  methanol  flow-rate 
is  written  as 

n 

Un  =  I<m  X  Y)fA(KpEh  I<dCEi )  X  Ts  (21 ) 

1  =  1 

where  Km  is  the  proportional  gain  and  I<p  and  I<d  are  new  tuning 
parameters.  Fig.  14(A)  and  (B)  shows  that  if  the  load  demand  is  ini¬ 
tially  changed  from  40  to  60  A,  then  both  stack  temperature  and 
the  hydrogen  pressure  can  be  regulated  at  the  desired  levels  by 
tuning  one  parameter  I<m  and  alternative  fuzzy  rule  bases,  respec¬ 
tively.  For  a  positive  75%  step  change  of  current,  Fig.  15  shows  that 
the  controller  configurations  are  employed  by  tuning  I<w  and  Km. 
Notably,  both  simulation  results  show  that  the  multi-loops  SISO 


Fig.  14.  MIMO  fuzzy  control  architecture  by  tuning  parameter  I<w  and  alternative  fuzzy  rule  bases:  (A)  stack  temperature:  (B)  hydrogen  pressure  at  the  anode;  (C)  corresponding 
manipulation  of  the  water  flow-rate;  (D)  corresponding  manipulation  of  the  methanol  flow-rate. 


928 


W.  Wu,  C.-C.  Pai  /  Journal  of  Power  Sources  194  (2009)  920-930 


Fig.  15.  MIMO  fuzzy  control  architecture  by  two  tuning  parameters  I<w  and  I<m  for  the  load  change:  (A)  stack  temperature;  (B)  hydrogen  pressure  at  the  anode;  (C)  corresponding 
manipulation  of  the  water  flow-rate;  (D)  corresponding  manipulation  of  the  methanol  flow-rate. 


control  design  guarantees  the  stable  output  regulation  of  a  heat- 
integrated  power  generation  system. 

5.  Conclusions 

In  this  work,  a  critical  conception  and  approach  from  the  process 
system  engineering  is  extended  to  construct  a  standalone  and  valu¬ 
able  power  generation  system.  The  heat  integration  is  based  on  the 
arrangement  of  HEXs,  the  effective  mixing  of  flowing  CH3OH  and 
02  at  the  inlet  and  countercurrent  water  flow.  Based  on  the  thermo¬ 
dynamic  properties  and  specifications  of  each  reactor  of  the  FPUs, 
the  inlet  and  outlet  temperatures  of  each  reactor  and  the  amount  of 


hydrogen  produced  at  the  exit  of  the  PROX  reactor  can  be  evaluated 
using  the  technical  computing  software  Matlab/Simulink™.  Fuzzy 
incremental  control  technique  within  the  multi-loop  feedback  con¬ 
trol  framework  can  be  applied  successfully  to  ensure  the  output 
regulation  performance  of  the  heat-integrated  power  system,  in 
spite  of  unknown  disturbances. 
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Appendix  A. 

Simulink™  blocks  for  PEMFC  stack  and  fuel  processing  subsystems  are  displayed  as  follows. 


AA.  A  grid  connected  PEMFC  stack  system 


A.2.  Heat-integrated  fuel  processing  units 


H20_ln 


Appendix  B. 

(i)  The  kinetics  of  methanol  reformer  (MR)  include: 
Steam  reforming 

Cat 

Me0H  +  H20^C02+3H2,  AH°  =  49.5kJ/mol 


I<eq,R  =  1-849  x  1010  exp 

£  i  p^ipk 

qR  Keq,RPMeOHPH20 

/  81000\ 

rR  =  6.75  exp  ( - ^MeOH^qK 


and  the  corresponding  reaction  rate  rR  (mol  g  cat-1  s-1)  is 
shown  by 


Partial  oxidation  of  methanol  (POM) 
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MeOH  +  l02-^C02  +  2H2,  AH°r  =  -192kJ/mol 

and  the  corresponding  reaction  rate  r0x,M  (mol  g  cat-1  s-1)  is 
shown  by 


,  n  /  65000 \ 

1<ox,m  =  0.466  exp  f - \ 

rOX,M  =  1<OX,MP  MeOH^o25 

Methanol  decomposition  (MD) 


Pm,  Cpm,  V,  £jout  =  Qj"  -  Q""  -  Qtr“ns/er 

dt  (B7) 

Pm2Cpm2V2  J-T°ut  =  -  or  +  4“/er 

where  notation  ‘m’  represents  the  mean  value.  Q}n  and  Qout 
are  the  inlet  and  outlet  flows  of  energy,  respectively.  The  heat 
capacity  of  the  component  is  listed  in  Table  4.  The  rate  of  heat 
transfer  (ftransfer  is  close  to  the  form  of 

Qtransfer  =  ECmin{Tin_Tm)  (B8) 


Cat 

MeOH  ^  CO  +  2  H2,  A  H°  =  90.61<J/mol 

and  the  corresponding  reaction  rate  rD  (mol  g  cat-1  s_1)  is 
shown  by 

Keq,o  =  1.718  x  1014  exp  (-^^) 

PcoPi 

(B3) 

/  76000 \ 

rD  =  1.12  exp  ( - J  ^MeOH^qD 

(ii)  The  kinetic  of  water  gas  shift  (WGS)  reactor  is  described  by 

Cat 

CO  +  H20^C02  +  2H2,  A  H°  =  -41 .15  kj/mol 

Referring  the  issue  [20],  the  equilibrium  constant  I<WGS  can 
be  expressed  as 

v  rWGS(hff^  +  r  WGS ) 

KWCS  =  (n^-rwcs)««0-rwcs) 

=  exp  (56B3,5  +  1 .077 ln(T)  + 5.44  x  10~4r 

-1.125  x  10“7r2  -  _  13. 148^  (B4) 

where  rWG s  (mol  g  cat-1  s-1)  represents  the  reaction  rate  of 
WGS  reactor. 

(iii)  The  kinetics  of  preferential  oxidation  (PROX)  reactor  include; 

CO  oxidation 

CO  +  102-^lc02,  A  H°  =  -282.99  kj/mol 

and  the  corresponding  reaction  rate  rCo  (mol  g  cat-1  s-1 )  is 
shown  by 

-rco  =  3.528  x  102  exp  (=^^)  P^co'1  (B5) 

H2  oxidation 

H2  +  ^02-^Ih20,  AH°  =  -285.84  kj/mol 

and  the  corresponding  reaction  rate  rH2  (mol  g  cat-1  s-1)  is 
shown  by 

-%  =  2.053  x  102  exp  (^§^)  P°0f  (B6) 

(iv)  The  model  of  heat  exchanger  (HEX)  is  considered  as  a  counter¬ 
flow  single-pass  heat  exchanger  [21  ],  the  outlet  temperature  of 
two  flows  is  modeled  as 


where  £  represents  the  heat-exchanger  effectiveness,  and  the 
minimum  heat  capacity  Cmin  is  expressed  as 


ril\CpmA1Tl2Cpm2 

hi  \Cprrii  +  hl2Cpm2 


(B9) 
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